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The aim of this work is to verify the new entropic and information inequalities for non-composite 
systems nsing experimental 5x5 density matrix of the qudit state, measured by the tomographic 
method in a multi-level superconducting circuit. These inequalities are well-known for bipartite and 
tripartite systems, but have never been tested for superconducting qudits. Entropic inequalities 
can also be used to evaluate the accuracy of experimental data and the value of mutual informa¬ 
tion, deduced from them, may charachterize correlations between different degrees of freedom in a 
noncomposite system. 
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INTRODUCTION 

During the last few decades tremendous progress has 
been made in experimental control over quantum sys¬ 
tems. In particular, experiments with superconducting 
circuits, based on Josephson junction devices [1, 2], have 
been rapidly developing recently [3]. Specifically spec- 
troscopical [4, 5] and time-domain [6] properties of such 
systems were studied both theoretically and experimen¬ 
tally. With the improvement of coherence time of su¬ 
perconducting qubits it became possible to obtain the 
density matrices of such systems, using quantum state 
tomography [7] as well as Wigner tomography [8]. 

Along with the development of quantum circuits, prop¬ 
erties of composite quantum systems, i.e. systems con¬ 
taining subsystems, have been extensively studied, which 
resulted in numerous practical applications. These sys¬ 
tems were also described in terms of classical information 
theory [9] in the quantum domain [10] and their informa¬ 
tion and entropic characteristics were investigated, in¬ 
cluding the von Neumann entropy and quantum mutual 
information, discord related measures, entropic inequal¬ 
ities, contextuality, causality, subadditivity and strong 
subadditivity conditions. 

On the contrary, the idea of using noncomposite quan¬ 
tum systems for quantum technologies was suggested [11- 
13] and quantum correlations in such systems have been 
analyzed only in recent times [14, 15]. The latter opened 
a way of mapping information and entropic measures for 
composite quantum systems on the noncomposite quan¬ 
tum systems [14-19]. 

In this work, we aim to verify the entropic and informa¬ 
tion inequalities using experimental 5x5 density matrix 


of the qudit state {j = 2), obtained using direct Wigner 
tomography in a superconducting circuit [8, 20, 21]. The 
inequalities were obtained using approach [14-18] to get 
analogs of subadditivity and strong subadditivity condi¬ 
tions, well-known for bipartite and tripartite systems, for 
a single qudit state. 

SUPERCONDUCTING CIRCUITS 

Superconducting circuits with Josephson junctions are 
macroscopic quantum objects, that can be several mi¬ 
crometers wide while still preserving quantum proper¬ 
ties. This happens because they are artficially isolated 
from the environment which leaves them with a single 
degree of freedom. The intrinsic parameters of these cir¬ 
cuits can be engineered as desired and adjusted with an 
external parameter (for example, a magnetic field), so 
they are thereby often called ” artihcial atoms”. 

Josephson junction 

The Josephson junction in superconducting circuits 
serves as a non-dissipative nonlinear element. It consists 
of two superconductors, separated by a thin insulating 
layer, through which Gooper-pairs can coherently tunnel. 
This system was described by Brian Josephson [22], who 
showed that supercurrent across the junction depends on 
the phase difference between the superconductors: 

/ = Jc sin ((/)2 - )= Jc sin (/), (1) 

where Iq stands for the maximum non-dissipative current 
flowing through the junction, i.e. the critical current. 
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Josephson also showed that when the voltage is applied 
across the junction the phase difference changes in time, 
which leads to the oscillations of the critical current with 
the angular frequency w. 

h^=huj = 2eV (2) 


When we substitute this into the time derivative of Eq. 
(1) and compare it to the Faraday’s law, we obtain the 
Josephson inductance: 




_ $0 

2e/c cos (j) 


( 3 ) 


As the Josephson junction has some intrinsic capac¬ 
ity C it behaves as a nonlinear oscillator with angular 
frequency ujp-. 
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( 4 ) 




FIG. 1: (color online) Tilted washboard potential (a) 
and quantized energy levels inside one of the potential 
wells (b). 



FIG. 2: (color online) The Josephson phase circuit 
(JPG) with an on-chip SQUID. The schematic diagram 
of the circuit. The left part corresponds to the JPG and 
the right part shows the on-chip SQUID, which is used 
for the readout, (b) The micrograph of the fabricated 
sample. (Images adapted from [20]) 


The physical implementation of this system is called the 
Josephson phase circuit [23, 24] and is shown in Fig. 2. 

The quantum state of the Josephson phase circuit is 
controlled via pulses of bias current. The measurement 
of the state employs the escape from the potential well via 
tunneling. For example, to measure the occupation prob¬ 
ability of state jl) one can pump microwaves at frequency 
W 41 , which will induce a jl) —>■ |4) transition. Then the 
state will rapidly tunnel due to the large tunneling rate 
r 4 . When the tunneling occurs, a voltage appears across 
the junction, which can be measured directly by an on- 
chip SQUID. 

In this paper we utilize the results, obtained in the 
experiment by Shalibo et al. [8, 20, 21], in which the 
Wigner distribution of the Josephson phase circuit was 
directly measured using simple tomography pulses. 


The total current flow trough the junction can be 
written as J = IcSincj) -\- V/R -\- CV. Substituting 
V = {h/2e)(j) from Eq. (2) we obtain: 

J =/csin())-t — —(5) 


which is equal to the equation of motion of a particle, 
moving in a tilted washboard potential: 


- 1 • 5C/(<^) ^ 


( 6 ) 


where U = — cos(j)^ 


is shown in Fig. la. 


Superconducting qudit 

A closer look at one of the wells in the tilted wash¬ 
board potential in Fig. 1(b) with the quantized energy 
levels gives us a perfectly suitable d-level system (qudit). 
Varying the potential by an external magnetic field, we 
can achieve a desired number of energy levels in the well. 


ENTROPIC INEQUALITIES 

Quantum states are generally described by the density 
matrix operator p, which has the following properties: 

Tr(p) = l, p = pt, p>o (7) 

We consider a 5 x 5 density matrix for a qudit (j = 2): 


f Pll 

Pl2 

Pl3 

Pl4 

Pl5^ 

P21 

P22 

P23 

P2A 

P25 

P31 

P32 

P33 

P3A 

P35 

Pil 

PA2 

P43 

Paa 

PA5 

\Psi 

Pb2 

P53 

P5A 

P55/ 


We can rewrite this as a 6 x 6 matrix, by adding one 
more zero row and zero column: 


/ 

Pll 

Pl2 

Pl3 

Pl4 

Pl5 

0 

\ 


P21 

P22 

P23 

P24 

P25 

0 



P31 

P32 

P33 

P3A 

P35 

0 



P41 

PA2 

P43 

Paa 

P45 

0 



P51 

P52 

P53 

P54 

P55 

0 


V 

0 

0 

0 

0 

0 

0 

/ 
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While looking at this system one can realize that it 
can be viewed as tensor product of two subsystems - a 
qubit and a qutrit. So, using an invertible mapping of 
indices 1 o -1 - 1/2; 2 -1 1/2; 3 O 0 - 1/2; 4 

0 1/2; 5 e-)- 1 — 1/2; 6 e-)- 1 1/2, we obtain the density 
matrix, which describes the bipartite qubit-qutrit state. 
The density matrices of the subsystems are generally de¬ 
rived by taking the partial trace over the corresponding 
indices. We propose a simplified approach by dividing 
the density matrix into blocks with fewer dimensions: 


/ 

Pll 

P 12 

Pl3 

Pl4 

Pl5 

0 ' 


P 21 

P 22 

P 23 

P24 

P 25 

0 


P31 

P32 

P33 

P34 

P35 

0 


P41 

P42 

P43 

P44 

P45 

0 


P51 

P52 

P53 

P54 

P55 

0 

V 

0 

0 

0 

0 

0 

0 


/ i?ll Ri2\ 
R 22 ) 


( 10 ) 


Then the density matrices of the subsystems are: 

_ /Tritii Tr i?i 2 \ _ /pll + p 22 + P33 Pl4-t P25'\ (' 11 ') 

^ yTr i?2i Tr R22J \ P 4 i + P52 P 44 + P 55 / 


P 2 = (.Rll + R 22 ) 


( Pll + P44 Pl2 + P45 Pl3\ 

P 2 I + P54 P 22 + P55 P 23 1 (12) 

P31 P32 P33 J 


Now we can take a look at correlations in our system. 
One of the most important correlation characteristics is 
entropy. In this work we deal with the von Neumann 
entropy [25]: 


5'Ar = —Trplnp (13) 


For the von Neumann entropy of the bipartite system 
one can write the subadditivity condition: 

Sp < S'pi -l- Sp2 

— Tr p In p < — Tr Pi In Pi — Tr p 2 In p 2 (14) 
and the mutual information equals: 


^bpl 


— Sp-i^ Sp2 Sp 


(15) 


Now we can repeat this process for another partition 
of the 6 x 6 density matrix: 


/ Pll 

P12 

Pl3 

Pl4 

Pis 

0 ' 

P21 

P22 

P23 

P24 

P25 

0 

P31 

P32 

P33 

P34 

P35 

0 

P41 

P42 

P43 

P44 

P45 

0 

Psi 

P52 

P53 

P54 

P55 

0 

V 0 

0 

0 

0 

0 

0 


( I’ll ri2 ri3 
r2i r22 r23 
r 3 i r 32 1-33 


(16) 


and get the density matrices of the subsystems: 

( Trrii Trri 2 TrriaX /pii + P 22 P 13 + P 24 Pi5\ 

Trr2i Trr22 Trr23 = P3i-H P42 P33 + P44 P35 
Trrgi Trr32 Trr33/ \ psi P53 P 35 / 


(17) 


P2 = (ril -t 1-22 I- r33) 


/ pll + P33 + P55 
V P21 + P43 


P12 + P34 
P22 + P44 


(18) 


So the subadditivity condition takes the form: 

Sp < S'pi -I- Sp2 

— Tr plnp < — Tr Pi Inpi — Tr p 2 lnp 2 (19) 

and the mutual information equals: 

hp2 = Sp, + Sp, - Sp (20) 

Next we add two more zero rows and columns to this 
matrix to get an 8 x 8 matrix. The system, described by 
this density matrix, can be divided into three subsystems 
(represented by 2 x 2 matrices) by the following mapping 
of indices: 


1 O - 1/2 - 1/2 - 1 / 2 ; 

3 O -1/2 1/2 - 1/2; 

5 O 1/2 - 1/2 - 1/2; 

7 o 1/2 1/2 - 1 / 2 ; 


2 o - 1/2 - 1/2 1 / 2 ; 
4 -1/2 1/2 1/2; 

6 O 1/2 - 1/2 1 / 2 ; 

8 O 1/2 1/2 1 / 2 . 


Here, we use the same approach of dividing the matrix 
into blocks to calculate the partial traces and get the 
matrices for the subsystems: 


/ 

0 

0 

0 

0 

0 

0 

0 

0 ' 


0 

Pll 

P 12 

Pl3 

0 

Pl4 

Pl5 

0 


0 

P 21 

P 22 

P 23 

0 

P 24 

P 25 

0 


0 

P31 

P32 

P33 

0 

P34 

P35 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

P41 

P42 

P43 

0 

P44 

P45 

0 


0 

P51 

P52 

P53 

0 

P54 

P55 

0 

V 

0 

0 

0 

0 

0 

0 

0 

0 

/ 

0 

0 

0 

0 

0 

0 

0 

0 


0 

Pll 

P 12 

Pl3 

0 

Pl4 

Pl5 

0 


0 

P 21 

P 22 

P 23 

0 

P 24 

P 25 

0 


0 

P31 

P32 

P33 

0 

P34 

P35 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

P41 

P42 

P43 

0 

P44 

P45 

0 


0 

P51 

P52 

P53 

0 

P54 

P55 

0 

V 

0 

0 

0 

0 

0 

0 

0 

0 


( 21 ) 


The density matrices that we are using hereinafter are 
the matrix of the second subsystem, i? 2 , and two joint 
matrices of the ’’qubit-qubit” subsystems, pi 2 and P 23 : 


P 12 


/Pll Pl3 Pl4 0 \ 

P31 P22 + P33 P3A p25 

Pil p43 P4A 0 ’ 

V 0 p52 0 p55/ 


( 22 ) 


P 23 = 


i?2 = 


fo 0 0 0 \ 

0 Pll + Pl 4 + P 41 + P44 P12 + P15 + P 42 + P45 P13 + P43 

0 P21 + P 24 + P 51 + P 54 P22 + P25 + P 52 + P 55 P 23 + P 53 

\0 P3I + P34 P32 + P 36 P33 

/ Pll + Pl4 + P41 + P44 PI3 + P43 

V P31 + P34 P22 + P33 + P25 + P82 + P55 


(23) 


(24) 


For this kind of tripartite system one can write the 
strong subadditivity condition [26]: 


Sp + Sn, < Sp„ + Sp,, 

-Trplnp- Tri? 2 lni ?2 < -Trpi 2 lnpi 2 - Trp23lnp23 (25) 
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VERIFYING EXPERIMENTAL DATA 

Next, we calculate the density matrices of the subsys¬ 
tems from the experimentally obtained 5x5 density ma¬ 
trix. This density matrix corresponds to the qudit, men¬ 
tioned in section Superconducting qudit, and was mea¬ 
sured in [8, 20, 21]. 


Using these matrices, we can calculate corresponding 
entropies and mutual information and test the subad¬ 
ditivity condition for different partitions in Eq.(14) and 
(19). Moreover, we can also change the position of the 
zero-row and zero-column in Eq.(lO) and (16) to see, how 
these entities will change. The results of these calcula¬ 
tions are given in the Table I and shown in the Fig. 5. 



FIG. 3: (color online) Experimentally obtained density 
matrix of a superconducting qudit [8, 20, 21]. 


The density matrices in Eq.(ll) and (12) are: 


Pi 


0.985 8.3-10-®-2.7-10-'‘i 

8.3-10-5 +2.7 0.006 


P2 


( 0.96 8.8-10-*-0.003 i 0.008-0.018 i 

8.8.10-* +0.003 i 0.004 -7.6 ■ 10-* - 2.9 ■ 10-* i 

V 0.008 + 0.0181 -7.6-10-* + 2.9-10-*i 0.026 


TABLE I: Calculated entropies and mutual information. 


Zero-row position 

5p 

Sbpl 

Shp2 

^bpl 

ibp2 

(i;i) 

0.1583 

0.300 

0.180 

0.1418 

0.0224 

(2; 2) 

0.1583 

0.1965 

0.3040 

0.0383 

0.1457 

(3; 3) 

0.1583 

0.1968 

0.3042 

0.0386 

0.1459 

(4; 4) 

0.1583 

0.2001 

0.1987 0.0418 

0.0404 

(5; 5) 

0.1583 

0.1873 

0.2059 

0.0291 

0.0477 

(6; 6) 

0.1583 

0.1996 

0.1768 

0.0413 

0.0185 



Position of the zero row 
(a) First partition 


Position of the zero row 
(b) Another partition 


FIG. 5: (color online) The entropies and mutual 
information from Table I, plotted versus the position of 
the zero-row, for the ” qubit-qutrit” partition. 


Analogously, for the Eq.(17) and (18) we obtain the 
following density matrices: 

/ 0.96 0.008 -O.OlSi -0.006 - 8.6 ■ 10-*i \ 

pi = \ 0.008 + O.OlSi 0.028 0.005 -0.007i , 

\ -0.006 + 8.6 ■ 10-* i 0.005 + 0.007 i 0.004 j 


Finally, we calculate the density matrices for the tri¬ 
partite system (Eq.(22), (23), (24)): 


Pl2 


( 0.959 0.008-0.0181 0.0002 -0.00041 

0.008 + 0.0181 0.026 0.003 + 0.00131 

0.0002 + 0.00041 0.003-0.00131 0.0018 

\ 0 - 0.0001 - 0.00021 0 


0 

- 0.0001 + 0.00021 
0 

0.004 




(0 

0 

0 

0 

\ 

/ 0.99 0.005- 0.002 i"! 

P23 = 

0 

0.961 i 

-0.005- 0.0041 

0.012-0.0191 


0.005+ 0.002 i 0.002 )’ 

0 

-0.005 + 0.0041 

0.004 

0.004 + 0.00641 




0.012 + 0.0191 

0.004-0.00641 

0.026 

/ 



(a) pi (b) p 2 



(c) h (d) P2 


FIG. 4: (color online) Calculated density matrices for 
the bipartite system: (a), (b) correspond to the first 
partition; (c), (d) correspond to another partition. 


R2 = 


0.961 0.012-0.019i 

0.012 +0.019 i 0.030 



(a) R 2 (b) pi 2 (c) P23 


FIG. 6: (color online) Calculated density matrices of 
the subsystems for the tripartite system. 

After calculations, the strong subadditivity condition 
(Eq.(25)) reads: 0.2997 < 0.3142; so the mutual informa¬ 
tion equals: I = Spj^^+Sp^^—Sp—Sn^ = 0.3142—0.2997 = 
0.0147. 
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CONCLUSIONS 


We have checked that the experimentally measured 
density matix of a superconducting qudit [8] satisfies 
the new entropic inequalities for non-composite systems, 
given by equations (14), (19) and (25). These inequal¬ 
ities can be further used to evaluate the accuracy of 
the experimental data. Moreover, the value of mutual 
information, deduced from entropic inequalities, may 
charachterize correlations between different degrees of 
freedom in a noncomposite system. There also exist 
other inequalities for the von-Neumann and q-entropy, 
which will be checked in future publications. 

We would like to thank E. Kiktenko and I. Khrapach 
for stimulating discussions. 
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